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investigated, and it is found that S n ,n-i is isometric to the hyperbolic 
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1 Introduction 



Many aspects of the dynamics of topological solitons of Bogomol'nyi type can be 
understood in terms of the geometry of the moduli space M^v of static iV-solitons. 
This approach, originally due to Manton, addresses such diverse issues as low 
energy soliton scattering, the thermodynamics of soliton gases, and the quantum 
mechanics of solitons. For a comprehensive review, see [7]. Mathematically, the 
main object of study is the L 2 metric 7, a Riemannian metric on which can be 
thought of as the restriction of the kinetic energy of the parent field theory. There 
are comparatively few situations in which explicit formulae for 7 are known, and 
one usually must make do with only partial or qualitative information. 

This paper considers one of the rare cases where explicit progress is possible, 
namely Ginzburg-Landau vortices moving on the hyperbolic plane. In this case, 
the Bogomol'nyi equations for static iV-vortices can be reduced to Liouville's 
equation on a disk, which is integrable. Strachan exploited this fact [TT] to 
obtain an implicit formula for 7 in terms of the analytic behaviour of the Higgs 
field near the vortex centres. From this he deduced explicit formulae for the 
metric on M] and M 2 , but the calculations become intractable for N > 3. In 
this paper, we will find exact formulae for the induced metric on certain totally 
geodesic submanifolds of M^r for all N, obtained by imposing invariance under 
certain symmetry groups. Each submanifold S n m C Mjy has (real) dimension 
2 and consists of static iV-vortex solutions wherein n single vortices occupy the 
vertices of a regular polygon, and < m < n — 1 coincident vortices sit at 
the polygon's centre (so m + n = N). The two dimensions correspond to the 
orientation and radius of the polygon. Geodesies in S n m are conjectured to 
correspond to low-energy iV-vortex scattering trajectories in the case of slow, 
rotationally equivariant initial data. 

We will discuss the curvature properties of S nm and show that S„ jn -i is 
isometric to the hyperbolic plane of curvature —3^;- This fact was already known 
(and is rather trivial) in the case n = 1, but its generalization to any n > 1 is new 
and extremely surprising. It follows that £„, jn -i is isometric to a hyperboloid of 
one sheet in (2 + l)-dimensional Minkowski space (one of the standard models of 
hyperbolic space). It turns out that, for m sufficiently close to but different from 
n — 1, S n m can still be isometrically embedded as a surface of revolution in IR^ 2 ' 1 ), 
and we construct the generating curves for some of these surfaces numerically. 
Geodesic motion on H njn _i can be understood very directly, and we obtain explicit 
formulae fixing the relationship between scattering angle and impact parameter 
for two-vortex scattering with a stationary third vortex at the origin (geodesic 
motion in E 2j i). Finally, we consider a variant of the moduli space dynamics 
recently proposed (for Euclidean vortices) by Collie and Tong [3] in which the 
vortices experience an effective magnetic field determined by the Ricci curvature 
of Mat. This flow is analyzed numerically on £2,0 an d exactly on E ni7l _i. 
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2 Vortices on the hyperbolic plane 

In this section we review Ginzburg-Landau vortices on the hyperbolic plane of 
curvature — \ [IHEj, which we denote H, with metric G. It is convenient to use 
the Poincare disk model of H, so H = {z G C : \z\ < 1} and 

G = Q dzdz, Q = ^— — . (2.1) 



Then (critically coupled) Ginzburg-Landau vortices on H are minimals of the 
potential energy 



V(4>, A ) = 7,f ( dA A * dA + d A0 A *d A + ^(1 - 00) 2 * U 



(2.2) 



where <fi : H C is a complex scalar field, A G fi 1 (H) is the gauge potential 
one-form, d^ = d<p — iA<fr, and * is the Hodge isomorphism. The standard Bo- 
gomol'nyi argument shows that, among fields satisfying the boundary condition 
|0| = 1 for \z\ = 1, with winding number N > 0, the potential energy V satisfies 

V > vriV, 

with equality if and only if 

{d A cf>){d/dz) = 0, (2.3) 

dA-*-(l-#) = 0. (2.4) 

Hence solutions of f)2.3p . f)2.4p minimize V in their homotopy class. Such solutions 
are called iV-vortices, and the zeros of are interpreted as individual vortex 
positions. 

Equations ( 12.31) and ( 12.41) can be reduced to a single gauge invariant equation 
by setting = e2 h+iX . One obtains 

4:—— + n-{le h = 0. (2.5) 

ozoz 

Setting h = 2g + 2 log |(1 — \z\ 2 ) the equation for h becomes Liouville's equation, 

4 H-^°- < 2 - 6 > 

which can be solved exactly. The solution is 

Sr=-l0gi(l-|/| 2 ) + il0 g 



df 



dz 
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(2.7) 
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where f(z) is an arbitrary, complex analytic function. With a simple choice of 
phase the scalar field is given by 

Then the first Bogomol'nyi equation (I2.3P is satisfied, if and only if 

A=-*dlog(i^t). (2.9) 

Note that <p vanishes at the zeros of 

To ensure that \<p\ — > 1 as \z\ — > 1, is nonsingular for \z\ < 1 and has 
winding number N, one must choose 

N+l / _ \ 

«*> - n s • (2-io) 

i=i x 7 

where ai,...,a/v+i are arbitrary complex constants with |aj| < 1. Naively, it 
seems that the moduli space of static iV-vortices should have complex dimension 
N+l, but this overcounts, since meromorphic functions f(z) and gauge equiv- 
alence classes of solutions of the Bogomol'nyi equations are not in one-to-one 
correspondence. In fact, the transformation 

(2.11) 

where \c\ < 1 leaves (</>, A) unchanged up to gauge. One can use this freedom to 
set a N+ i = 0, so 

N / - \ 

'M-n(SS), (2 ' 12) 

where \a$\ < 1. Note that the denominator of this rational map is uniquely 
determined by its numerator. So, hyperbolic iV-vortices are in one-to-one corre- 
spondence with degree N polynomials 

N 

P(z) = l[(z-a l ) (2.13) 

i=X 

all of whose roots lie in the open unit disk, and one concludes that M n has real 
dimension 2N. The coefficients a, are not natural complex coordinates on M^, 
however, as we shall see. 

It is useful to think of / as a holomorphic mapping of the Riemann sphere CU 
{oo} to itself. Clearly / has exactly 2N critical points, counted with multiplicity. 
Note that / commutes with the involution t : z \— > z~ l (reflexion in the equator) 
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so if z is a critical point of /, so is l(z). Since / has no critical points on 
the equator, \z\ = 1, it follows that it has exactly N critical points inside the 
unit disk, and iV outside. Let us denote the N (not necessarily distinct) critical 
points inside the unit disk Z r , r — 1, 2, . . . , N. These are interpreted as the vortex 
positions. They provide local complex coordinates on Mat\Atv, where An is the 
coincidence set, that is, the set of N- vortices for which the vortex positions are not 
all distinct. Good global complex coordinates are provided by the coefficients of 
the polynomial Ylf = i{z—Z r ). This equips M^v with a canonical complex structure 
which coincides, off Ajv, with the structure defined by the coordinates Z r . Note 
that the critical points Z r of / depend non-holomorphically on the parameters 
a;, so this canonical complex structure is different from the complex structure 
defined by the coordinates a» on M^. 

There is a natural Riemannian metric 7 on Mjv defined by restricting the 
kinetic energy of the abelian Higgs model 



to TMjy. We here choose to work in temporal gauge (A = 0), so must impose 
Gauss's law 



as a constraint on (<9o0, 8qA), where 5 = — * d* is the coderivative on (H, G). It 
is known that 7 is Kahler with respect to the canonical complex structure just 
defined. This follows from a formula for 7 first obtained by Strachan [11] and later 
generalized and reinterpreted by Samols [10] and Romao [8]. The formula gives 7 
on Mjv\Aiv. In the next section, we will need to compute the induced metric on 
certain totally geodesic submanifolds E njm C M n+m which (for 2 < m < n — 1) lie 
entirely in A n+m , so we must make a small modification to Samols' calculation. 

Rather than compute 7 on Mat\Aat, we fix the positions of m of the vortices 
(not necessarily distinct) and allow only the other n = N — m vortices to move, 
along trajectories Zi(t), . . . , Z r (t) which remain distinct from one another, and 
the fixed zeros. This, then, defines an n-dimensional complex submanifold of 
M n+m , which we denote by M$ z \ where p(z) is the unique monic degree m 
polynomial whose roots are the fixed zeros (so = M n ). Samols' argument 
generalizes immediately to this setting, and one finds that the kinetic energy of 
a trajectory in is 



Here b r are coefficients in the expansion of h = log \ <ft\ 2 around the zero Z r of </>, 




(2.14) 



Sd A = -(<f)d <j) - 0<9 o </>) 



(2.15) 




(2.16) 




5 



On hyperbolic space, |0| 2 is known explicitly, so the dependence of b r on the vortex 
positions Z 1; . . . , Z n can also be determined explicitly, in principle. In practice, 
as we shall see, this is very difficult, but some properties of the coefficients b r are 
immediate. First, since T is manifestly real, one must have 

(2.18) 

dz s dz r v ; 

It follows that the metric on induced by T, 

n { db \ 

7 = 7i \ V{Z r )6 rs + 2^ dZ r dZ s , (2.19) 

r,s=l ^ r ' 

is Hermitian with respect to the canonical complex structure, with Kahler form 

• n / db \ 

uo = -7T V ( Q(Z r )5 rs + 2— f ) dZ r A dZ s . (2.20) 

2 ^— ' V "^r / 

r,s=l x ' 

Clearly, du; = 0, by (12.181) . that is, 7 is Kahler. Following Romao, it is convenient 
to define the (0, 1) form 

n 

b = ^b r dZ T (2.21) 

on M„ 2 ', which is known to satisfy db = 0, [8]. The Kahler form on MS' 2 ' may 
then be compactly written 



1 

UJ = —TV 

2 



2db + tt(Z r )dZ r A dZ r j . (2.22) 



r=l 



3 Totally geodesic submanifolds of M 



N 



We shall use formula (I2.22j) to deduce the induced metric on certain totally 
geodesic submanifolds of which we now define. Recall that points in Mat 
are in one-to-one correspondence with monic degree iV polynomials with roots 
only in the unit disk. We can choose to identify a 7V-vortex with the polynomial 
whose roots are the vortex positions, but it is somewhat more convenient to iden- 
tify it instead with the monic polynomial P{z) defined by the numerator of the 
rational map f(z), whose coefficients we denote $ G C, 

P( z ) = z K + p lZ "-i + ... + f] N . (3.1) 

There is an obvious action of U(l) on such polynomials, whose associated action 
on Mjy is clearly isometric (it just rotates the vortex positions about z = 0), 
namely, 

P{z) 1 — * \- N P(\z) (3.2) 
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for each A G U(l). In terms of the coefficients of P(z), the action is 



{Pi, 132,..., (3 N ) » {\~ x (3i, \~ 2 p 2 , . . . , \~ n Pn). (3.3) 

From now on we consider the case where A = e 27U//n , for n some positive integer, 
so that A generates the cyclic group C n . Then P{z) is a fixed point of A if and 
only if flk — for all k not divisible by n. So the fixed point set of A consists of 
all polynomials of the form 

P(z) = z N + (3 n z N ~ n + (3 2n z N - 2n + ■■■ + (3 pn z N - pn , (3.4) 

where p = | — ] . The corresponding submanifold of Mjy has complex dimension 
p and is totally geodesic by a standard symmetry argument In particular, if 
Y < n < N, so that p = 1, the fixed point set has complex dimension one and 
consists of polynomials 

P(z) = z m (z n + f3 n ), \(3 n \<l, 0<m<n. (3.5) 

Let us denote this totally geodesic submanifold S n m C M n+m . In order to cal- 
culate the metric on S n m we first need to find the vortex positions, that is, the 
critical points of 

which arises from (13. 5p (replacing (3 n = —a n for later convenience). Recall / has 
exactly n + m critical points, counted with multiplicity, in the unit disk H. Now 
f(Xz) = X m+l f(z), so z is a critical point of / if and only if Xz is a critical point. 
Hence, / has n critical points Z r at the vertices of some regular n-gon, 

Z r = a X r ~\ (3.7) 

and the other m critical points must be coincident at (the only fixed point in 
H of A). So S nm consists of vortex configurations with m vortices coincident at 
2 = and n vortices located at the vertices of a regular polygon centred on 0. 

We next determine how the vortex position a = Z\ G C is related to the 
complex parameter a. If a is real positive, then / clearly has a real postive 
critical point, and / remains unchanged under a i— > Xa. It follows that 

a = av(\a\) (3.8) 

where v is a positive real function of \a\ only. A short calculation leads to a 
quadratic equation in v n . We choose the solution which satisfies \a\ < 1 provided 
\a\ < 1 and obtain 



(m + 1) (1 + |d 2n ) - \ (m + l) 2 (1 - \a\ 2n Y + An 2 \a\ 2n 

v n = — I , ' — (3.9) 

2a 2n (m + l-n) v ; 
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for n > m + 1. For n = m + 1 the solution simplifies considerably, namely. 

2 



2n ' 



(3.10) 



We can now eliminate the parameter a from f(z) in favour of a which, being 
one of the vortex positions, is a good (local) complex coordinate on S n m with 
respect to the canonical complex structure: 

/(*) = - , { n - n ■ (3-H) 
1 — v n a n z n 

By the Kahler property and rotational invariance, we know that the metric on 
S n m must take the form 

7 = F(\a\)dada (3.12) 

for some conformal factor F(|a|). To compute F, we require information on the 
coefficients b r , and hence a formula for h = log \<fr\ 2 . Now f'(z) has zeros at X r a 
and A r a _1 (recall critical points of / are paired z «-> z -1 ), and a zero of order m 
at 0. It follows that 

s _ (m + l)^"z m (z" - q")(l - a n z") 

* _ 7l „-„ „\9 • (3.13) 

(1 — v n a n z n y 

In order to avoid the logarithmic singularities of h near z — Z\ — a, we define 
the regularized version of (I2.17P as 

h reg = log |0| 2 - log (2 - a) - log (z - a) . (3.14) 

Since a is a simple zero we can calculate the coefficient 61 in (12.171) as 



bx 



2dh reg 



dz 

which leads to 



(3.15) 



1 / 1 2na n a n 4nu n a n a n 4aa \ 

6i = - 2m + n-l-- — + - — -z • 3.16 

a \ 1 — a n a n 1 — u n a n a n 1 — aa J 

We have calculated b\ only, but we can deduce b r for r > 2 by rotational 
symmetry. From the definition of the coefficient b r in f)2.17p . it follows that if we 
simultaneously rotate all the zeros by A' G U(l), the b r coefficients transform as 

6 r (A%, . . . , X'Z n ) = X'b^, ...,Z n ). (3.17) 

In the case A' = A = e 2nl ^ n , this rotation just cyclically permutes the vortices, so 
we see that 

b r (a) = X r ~%(a). (3.18) 
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It also follows from (13.171) that ab±(a) is a function of |a| only, which is a consis- 
tency check on our formula (13.161) . 

To complete the calculation of the metric, we note that E njm (with the point 
a = removed, strictly speaking) is a submanifold of M^ m , on which the formula 
( 12.221) for the Kahler form holds. Let us introduce the real vector fields^ 

r=l x / r=l x y 

on M^ m , and note that on S„,. m these coincide with 



X=\a\ J — Y=—, (3.20) 
o\a\ dip 

where a = \a\e 1 ^. Hence, the conformal factor we seek is 

F(\a\) = -^ 2 u(X,Y) = ^Udb)(X,Y) + n\a\ 2 n(\a\)) 
\a\ z \a\ z \ J 

= j^ 2 (x[b(Y)]-Y[(b(X)])+n-Kn(\a\) (3.21) 
since [X, Y] = 0. Now 



b(X) = b r dZ r (X) = Y Z r b r (3.22) 



r=l r=l 



and b(Y) = b(JX) = —ib(X) since b is a (0, 1) form. Since X and Y are 
tangential to £ n , m , to compute = —iX[b(X)] and Y[6(X)] it suffices to 

know b(X) only on S„ im , where Z r = \ r ~ l a. But using ( 13. 18ft we see that, on 
£ m , n , b(X) = nabi(a) which, as we have remarked, is a function of |a| only. 
Hence, F[6(X)] = 0, and we find that 

F(\a\) = Tin (n(\a\) + -L^—^hia))) . (3.23) 
\ \a\ a\a\ ) 

The derivation for this formula used only rotational invariance and ( 12.221) . so 
it holds equally well for vortex polygons on the Euclidean plane, or any other 
surface of revolution. 

Now, we can evaluate the metric using equations (I3.23|) . (13.161) . (13.91) and 
( I3.10p . obtaining 



Ann 3 ] 


n 


2n-2 


(1- 


o 


\2ny 



2n(l + | 


a 


2n-j 


y/(m + 1) 2 (1- 


a 


|2n)2 + 4n 2 


a 


In 



7 = „ ' , 1 + n . == i = . ,. nE dada (3.24) 



1 In [2J, Chen and Manton used the vector field X to derive an interesting integral formula 
for the Kahler potential. 
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for n > m + 1. For n = m + 1, the metric simplifies and we have 

127m 3 |a| 2n - 2 , ,_ . 
7 = "Ti 1 i9„xo aacta. (3.25) 



In the case m — 0, equation ( 13. 24ft agrees with previous results of Strachan [TT] 
for n = 1,2, and Krusch and Sutcliffe [5] for general n. For |a| close to 1, the 
moving vortices are very far apart, so one expects the metric to approach the 
metric induced on £ m ,n by the product metric on (Mi) n 

n yiix 

7product = ^ (I _ \Z \2)2 dZrdZr > (3.26) 
r=l I r l / 

which is 

127T7 

_ 12^2 



127m 

7oo = 7: 7—r^-rdada, (3.27) 

1 — \a\ 



the metric on the hyperbolic plane of curvature n = Indeed, the full metric 

is asymptotic to 7^: 

/ 127rn \ 
7 = ( 7 i — ToTo — nm(m + 2)ir + 0(1 — lad) ) dada. (3.28) 

\(i-|«lr / 

Much more surprising is the fact that is exactly isometric to the hyperbolic 

plane of curvature k = —3^. To see this, one should introduce the global complex 
coordinate x = 0/11 1 which is in one-to-one correspondence with points in E n m 
(recall that a 1— > Xa just permutes the vortex positions, so maps each point in 
£n,m to itself). The metric is then 

47m / 2n (1 + I yI 2 ) 

7 = - , |2 . 2 1 + ! 1 1X1 ; I d X dx- (3.29) 

tl-IXl; \ ^(m + l) 2 (1 - |x| 2 ) + 4n 2 |x| 2 

Note that this is nondegenerate at x — 0. For m = n — 1, it simplifies to 

12nn dxdx 

7 = n , * 2 , (3-30) 

which, as promised, is the metric on the hyperbolic plane with curvature k = 
i_ 

3irn ' 



4 The geometry of S n m 

It follows immediately from (I3.29P that S n m is geodesically complete and has 
infinite volume. The radial curves x{t) = t\, X E U(l), — 1 < t < 1 are pre- 
geodesics corresponding to the dual polygon scattering trajectories characteristic 
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of topological solitons in two dimensions. It is straightforward to compute the 
Gauss curvature of the surfaces S m n using the formula 



2\x\F(\x\)d\x\ \F(\x\)d\x 

In all cases 

lim k = — - 



IxHl 37T71 



and for m close to n — 1 , k is uniformly negative. As already noted, for m = n — 1 , 
k = — ^ for all x- By c 
example, on S n one has 



k = for all x- By contrast, for m small, k is positive close to x — 0- For 



Tin (1 + 4n + 4?t/) 

which is positive for all n > 2. 

For purposes of visualization, it would be useful to find an isometric embed- 
ding of S m>n as a surface of revolution in M 3 , along the lines of the "rounded 
cone" picture of the Euclidean two- vortex moduli space discussed by Samols [TO] . 
In our case, each £„ im is asymptotic to a complete space of constant negative 
curvature, so it is clear that no such isometric embedding exists. However, we 
certainly can find an isometric embedding of S n n _! as a surface of revolution in 
l^ 2 ' 1 ) (that is, M 3 equipped with the Lorentzian metric dx\ + dx\ — dx"^), invariant 
under rotations about the timelike x% axis. The image of this embedding is the 
hyperboloid on which 

x\ + x\ — x\ = Snn, x 3 > 0. (4.3) 

Of course, this is just the hyperboloid model of the hyperbolic plane. One expects 
to be able to generalize this to S nj?ri at least for m close to n — 1. A short 
calculation shows that any surface of revolution which intersects the symmetry 
axis must have k < at the intersection point, so isometric embeddings of E n)m 
for m small certainly do not exist, as k(0) > and x = is a fixed point of the 
rotation action. Using a straightforward modification of the method described in 
[6], one can obtain generating curves for the embedded surfaces of revolution for 
m close to n — 1 . Some examples are depicted in figure [TJ 

The geodesic flow on E 2 ,o was considered in detail in [11] . It is a simple matter 
to analyze geodesic motion in E ni7l _i for any n > 1, since each of these spaces is 
isometric to a hyperbolic plane, and geodesic flow is invariant under homothety. 
Hence, the geodesic trajectories in £ n , jn _i are the standard geodesies in the hyper- 
bolic plane. In the Poincare disk model, these are circular (or straight) arcs which 
intersect the boundary of the unit disk orthogonally. The vortex trajectories in 
physical space H corresponding to a geodesic in £ nn _i are the preimage of the 
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Figure 1: Generating curves for the spaces £ 8 , m iso metrically embedded as sur- 
faces of revolution in IR( 2,1 ) for m = 5, 6, 7 (solid curves, top to bottom). The full 
surfaces are swept out by rotating about the x 3 axis. The m = 7 curve coincides 
with the hyperbola x\ — x\ = —24n. The dashed lines mark the light cone. No 
isometric embeddings of S 8 m for m = 0, . . . , 4 exist since these have positive 
Gauss curvature at the fixed point of the rotation action. 

geodesic under the map z \— > z n of the unit disk to itself. Since this map is confor- 
mal, the vortex trajectories also intersect the unit circle orthogonally. Consider 
the case n = 2 in detail. Without loss of generality, we may restrict attention 
to the geodesies in T,2,i which intersect the boundary at e ±l ^' where < ip < ~. 
Each geodesic x{t) corresponds to a 3- vortex motion in which one vortex remains 
stationary at the origin and the other two move towards each other, scatter and 
recede, following the trajectories z(t) = ±x(t)a. We may derive an explicit for- 
mula for the scattering angle and impact parameter b associated with the i/j 
geodesic as follows. First we determine the geodesic in H which makes second or- 
der contact with the incoming vortex trajectory where it intersects the boundary, 
at e l ^l 2 . Simple trigonometry shows that this is a circular arc of radius ^tan^. 
This is the path that this vortex would follow if it did not interact with the other 
two. Call its exit point e l *. Then the scattering angle is the angle one must rotate 
the disk so that this free exit point shifts to the actual vortex exit point e - "^ 2 , 
that is, such that £ + = —%. One can similarly construct the initial free path 
of the other moving vortex: it is the image of the first under a rotation by tt about 
the origin. We define the impact parameter to be the hyperbolic distance (with 
respect to metric G on H) between these two free trajectories. The geometry is 
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Figure 2: The geometry of vortex scattering: the solid curves are the vortex 
trajectories in physical space, while the dashed curves are the free trajectories 
which the vortices would describe if they did not interact with one another. By 
definition, these are the geodesic arcs which have second order contact with the 
incoming trajectories at their entry points on the boundary of H. 



summarized in figure [2j Straightforward calculation shows that 

m = 4V2tanh- 1 ( ^ = Shl ^ 

\ 2 cos ip 

e(if>) = 2 tan" 1 ( -tanV J - V>- ( 4 - 4 ) 



Note that tjj — | gives the expected 90° head-on scattering process (i.e. b = 
and O = |), and that decreases monotonically towards as b increases. 



5 The Collie-Tong flow 

Motivated by a supersymmetric extension of the Abelian Higgs model with a 
Chern-Simons term, Collie and Tong have derived a new geometrically natural 
moduli space dynamics for Euclidean vortices [3]. This geometric flow is well 
defined on any Kahler soliton moduli space, and it is interesting to consider 
the dynamics of hyperbolic vortices under this flow, in comparison with their 
Euclidean counterparts. Recall that on a Kahler manifold M one has a closed 
2-form p constructed from the Ricci curvature R and almost complex structure 
J in the same way that the Kahler form is constructed from the metric and J, 
that is, 

p(X,Y)=R(JX,Y). (5.1) 
A trajectory x(t) G M is a solution of the flow if 

V d/dt x - = (5.2) 
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where A is a real constant, j} denotes the metric isomorphism T*M — * TM and i 
denotes interior product, (txp)(Y) = p(X,Y). A very similar flow has, in fact, 
been derived previously by Kim and Lee [4], but it is the geometrically natural 
formulation given above which is significant for our purposes, and this is due to 
Collie and Tong. For this reason, we will call equation (15.21) the Collie- Tong flow. 
It follows from (I5.2p that 



so this flow, like geodesic flow, conserves speed \\x\\- Given a solution x{t) of 
( 15. 2p . x(^ot) satisfies ( 15. 21) with parameter A' = A/Ao, so we can rescale A to 
any convenient value. Note that, in contrast to geodesic motion, the trajectories 
depend on ||x(0)||, not just the direction of x(0). If A = 0, one recovers geodesic 
flow, so one expects trajectories under (15.21) to approach geodesies when the initial 
speed is taken to infinity. 

Collie and Tong investigated the qualitative properties of this flow on the 
moduli space of centred Euclidean 2- vortices [3], finding scattering trajectories, 
bound orbits and bound orbits with Larmor precession. In this section we will 
make an analogous analysis of the flow on E 2i o numerically, and exactly. 
First, we note that £ n)m is a totally geodesic complex submanifold of M n+m , so 
its Ricci form coincides with the restriction of the Ricci form of M n+m to TS n m . 
Hence, solutions of the flow on E„ )m are solutions of the flow on M n+m . Now on 
any two-dimensional Kahler manifold, p =\klo where Kj IS, clS before, the scalar 
curvature, and u is the Kahler form. Hence (15.21) (with A = 2) becomes 



We have solved this equation numerically on S 2j o for various initial data using 
the ODE solver package in Matlab. The corresponding vortex trajectories are 
depicted in figure [3j Appropriate choices of initial data yield scattering trajecto- 
ries, bound orbits and Larmor precession, as predicted by Collie and Tong in the 
Euclidean context. One can also find complicated cycloid-like trajectories. 

We can make a more complete analysis of the flow on S„ ;re _x, since this is 
isometric to the hyperbolic plane of curvature — l/37m, so k is constant. By 
rescaling t we can scale away n in ( 15 .4p . so it suffices to consider the flow on, say, 
the hyperbolic plane of curvature —1. In the upper half plane model this has 
metric 



d_ 

It 



(X, X) = 2A(x, UxP) = 2X p(x, X) = 0, 



(5.3) 



^d/dtX - k,Jx = 0. 



(5.4) 



7 = — (dx 2 + dy 2 ). 

y 

To translate back to the disk model, we note that 



(5.5) 



z — i 



(5.6) 



X = 



-iz + 1 



z = x + iy. 
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(b) 



(d) 






Figure 3: Vortex trajectories under the Collie- Tong flow on E 2) o- The curves show 
the trajectory of one of the vortices in each two vortex motion, and illustrate a 
variety of behaviours: (a) scattering, (b) an almost-circular bound orbit, (c) 
Larmor precession, (d) a cycloid-like bound orbit. 

The vortex trajectories are then given by the n-th roots of x{t)- I n the z coordi- 
nate system (15. 4p becomes 



x xy + y 

y 



y + -{x 
y 

which has two conserved charges, 

1 



y ) -x 



E 



P 



0. 
0. 

i-y 
y 2 



(5.7) 



(5.1 



Consider the trajectory with initial data z(0) = i, z(0) = iv for v > 0. This has 
P = — 1, so its projection to the (y, y) phase plane is the energy level curve 



yV-(i-y) 2 )- 



(5.9) 



For < v < 1 the trajectories are periodic, with 1 — v < y(t) < 1 + v while 
for v > 1 the trajectories are unbounded, escaping from and to the boundary at 
infinity (y = 0) as t — > ±oo. Between critical points of y(t), we can determine 
the trajectory x(y) by solving 



dx 
dy 



x 

y 



l 



yf 



(5.10) 



One finds that 

(x + vf + (y- If = v 2 (5.11) 

so the trajectory is a circle of radius v centred on (—v, 1), or, if v > 1, the portion 
of this circle in the upper half plane. Note, in the latter case, that the trajectory 
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(b) (c) (d) 




Figure 4: Vortex trajectories under the Collie- Tong flow on S 21 . One vortex 
moves along the trajectory a(t) depicted, another moves along —a(t) while a 
third vortex remains stationary at the origin. In each case a(t) 2 describes a circle 
with centre in the unit disk. In case (b) this circle passes through so the vortices 
undergo 90° scattering at the origin. The circle captures in cases (c) and (d), 
but not in case (a). Note that all trajectories are periodic or escape to infinity. 

intersects the boundary acutely, not orthogonally, though the trajectories tend 
to a geodesic as v — > oo, as one would expect. By acting on this one-parameter 
family of circles with the isometry group of the hyperbolic plane (SL(2, M.) acting 
by fractional linear transformations, in the upper half plane model), we obtain 
all circles centred in the upper half plane (note that geodesies have centres on the 
boundary y — 0). Under the Mobius transformation (15.61) . these are mapped to 
circles with centres in the unit disk, and every such (arc of a) circle is a trajectory 
of the flow. The corresponding vortex trajectories in H (in the disk model) are 
then obtained by taking n-th roots. Clearly the behaviour is much simpler than 
that found numerically in £2,0 : all trajectories are either periodic or unbounded 
(both as t — > oo and t — > — oo). Some examples are depicted in figure HI 

It would be interesting to make a systematic analysis of the flow (15.21) on a 
general Kahler manifold. One expects, for example, that the flow is complete if 
and only if the manifold is metrically complete. A natural generalization of the 
case S n n _! considered in detail above would be to assume that M is a homoge- 
neous Einstein manifold, for example M = CP N . This is potentially relevant to 
vortex dynamics on the two-sphere close to the Bradlow limit [I]. 
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